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Abstract 

We address the problem of the optimal quantum estimation of the coupling param- 
eter of a bilinear interaction, such as the transmittivity of a beam splitter or the 
internal phase-shift of an interferometer. The optimal measurement scheme exhibits 
Heisenberg scaling of the measurement precision versus the total energy. 



1 Introduction 



How effectively may we estimate the strengtli of a simple interaction sucli as 
tlie transmittivity of a beam splitter or the phase-shift imposed in the internal 
arms of an interferometer ? The Hamiltonian describing the bilinear coupling 
between two bosonic modes has the form 

H = K{a)h + h'^a) , (1) 



where k depends on the specific interaction under consideration. By using the 
Schwinger representation of the SU{2) Lie algebra, with generators 

J. = ^{a^b + b^a) , Jy = l-{a^b-b^a) , J, = ^{a^a - b^b) , (2) 



we can rewrite the Hamiltonian as H = 2k,Jx and the evolution operator as 

= exp(-iJ^V') , (3) 



where the global coupling constant ip is equal to 2KAt, At being the effective 
interaction time. The evolution in Eq. (3) describes, for example, the interac- 
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tion of two light modes in a beam splitter with transmittivity r = cos^ ip [1] , or, 
apart from a fixed rotation, the evolution of the arm modes in a Mach-Zehnder 
interferometer, with ip representing the phase-shift between arms [2,3]. If the 
initial preparation of two modes is described by the density matrix po the 
evolved state in the interaction picture is given by 



In this paper we devote our attention to the estimation of ip through measure- 
ments performed on p^. We denote the generic POVM for such estimation 
process by dt/(0), such that the outcomes of the measurement are distributed 
according to 



where p{(f>\ip) represents the conditional probability of registering the outcome 
(f) when the true value of the parameter is ip. 

Our objective is to find the best strategy, i.e. the POVM that provides the 
optimal estimation of the parameter ip [5]. Since ip is manifestly a phase-shift, 
we can use general results from phase estimation theory, which provides the 
optimal POVM to estimate the phase-shift induced by a phase generator, i.e. 
a selfadjoint operator with discrete, equally spaced, spectrum. The optimality 
criterion is given in terms of the minimization of the mean value of a cost func- 
tion [5] that assesses the quality of the strategy, i.e. it weights the errors in 
the estimates. Since a phase-shift is a 27r-periodic parameter, the cost function 
must be a 27r-periodic even function of {ip — (f)), i.e. it has Fourier expansion 
with cosines only. The appropriate concavity properties of the cost function 
are reflected by expansion coefficients which are all negative apart the irrele- 
vant additive constant. A cost function with such Fourier expansion has been 
firstly considered by Holevo [4] , and for such reason it is usually referred to as 
belonging to Holevo's class. Notice that the optimal POVM for a given state 
Po is the same for every cost function in this class. A more general quantum 
estimation approach for different kinds of phase-shift and general quantum 
system is given in Ref. [6]. There, it is also shown that the kind of problem 
we are presently dealing with is, in a certain sense, the best situation, since 
the spectrum of our phase generator is the whole set of integers Z, including 
negative ones. In this case, there is an optimal orthogonal projective POVM, 
which can be regarded as the spectral resolution of a self-adjoint phase opera- 
tor. However, if the estimation is performed with the constraint of bounded or 
fixed energy, the optimal POVM and the optimal input state po do not corre- 
spond to a canonical quantum observable scheme [7]. Moreover, in general the 
optimal POVM depends on the input state, and the optimal POVM of Ref. 
[4,5] holds only for pure states, whereas a generalization to a class of mixed 



p^ = exp(-iJa;V')poexp(iJa;V') . 



(4) 



(5) 
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states, the so-called phase-pure states, has been considered in Ref. [6,8]. 

The optimal POVM, in the sense described above, provides an unbiased esti- 
mation of ip for preparation po of the two modes, in formula 

{(p) = ip with {(j)) = j (j) Tr[p^ du{(j))] Vp^ = U^PoUl . 



On the other hand, wc want also to find the optimal state po for the estima- 
tion of ip according to some cost function, which quantifies the noise of the 
estimation. The customary root mean square is not a good choice for a cost 
function, since the function (0 — ■0)^ is not 27r-periodic. A good definition for 
the precision of the measurement is given by the average of the cost function 
C{(j)—ip) = 4 sin^(^^), i. e. a "periodicizcd variance" , which obviously belongs 
to the Holevo's class. If the estimates occur within a small interval around the 
true value of the parameter ip, one has approximately C{(1) — tp) ~ (0 — tpY, 
whence 5'^ — \fC can be assumed as a reasonable measure of the precision of 
the measurement. 



2 Optimal estimation of the coupling parameter 

In order to solve our estimation problem, let us consider the following unitary 
transformation 

U = exp | — — (a^6 — 6^a)| = exp j^i'^-^j 




Using Eq. (6) we may rewrite Eq. (5) in the more familiar form of rotation 
along the z-axis 



V 



Tr \Up^U^ Wdz/(0)W^ 



Tr 



exp{—iJzip)UpoU^ exp{iJzilj) Wdz/(0)W^j 



: Tr [exp{-U^ip)Ro exp(i J^V') dp(0)] , 



(7) 



where we used the identity UJxU^ — Jz- Equation (7) shows that the problem 
of estimating the shift generated by on the state po is equivalent to that of 

estimating the same shift generated by on the rotated state Rq = UpqW . 
In particular, any POVM dz/(0) to estimate the J^-induced shift can be writ- 
ten as dvi^cj)) = Wdp{(j))U, where dp(0) is a POVM for the J^-induced shift 
estimation. 
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For pure states Rq = \i/-^o)) {{ipo\ (in the following we use double brackets 
for two-mode vectors) the degeneration of the spectrum of a'^a — b% can be 
treated using the technique introduced in [6] , and the optimal POVM for cost 
functions in Holevo's class [4] is proved to be of the form 



M4>) = ^ \EM{EA 



(8) 



with the vectors \E^)) given by 



(9) 



The vectors \d)) are certain eigenvectors oi D = a^a — b^b built by picking up, 
in every eigenspace Ti^ of the eigenvalue d, the normalized vector parallel to 
the projection of \ipo)) on Tid- In order to be more specific, let us consider an 
input state of the form 

oo oo oo 

m = E 0)) + E E K'V, d)) + i^i-'^\n, -d))] , (10) 



n=0 



n=0 d=l 



where \n,d)) is given by 



\n, d)) = < 



\n + d)a\n)b if d > 
|n)a|n — d)b if (i < 

The projection of |'0o)) in "^d is equal to 



n=0 



fill 



(12) 



such that the eigenvector \d)) reads as follows 

J:^,clf^\n.d)) 



id) 



(13) 



and the input state can be rewritten as 



lV'o)) = E7dM)) 7d 



\ n=0 
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Notice that the dependence of the POVM on the state j^'o)) is contained in 
the vectors \E^)). 

By adopting C{4> — i/j) as a cost function the average cost of the strategy 
corresponds to the expectation value of the cost operator C = 2 — E+ — E_, 
where the raising and lowering operators and E_ are given by 

E^ = J2\d+l)){{d\ E_ = Ei 
(with the vectors \d)) defined as above). 

The optimization problem is that of minimizing the average cost of the strategy 
C^l^ / ^ C(</. - V^) = Tr [Ro C] ^ {{^o\2 - E+ - E.\^o)) 



(14) 

with the constraint that the solution is a normalized state. The Lagrange 
function is given by 

C ^ C - XiiM^l^o)) , (15) 

with A being the Lagrange multiplier for the normalization constraint. The 
solution of this problem is a state with infinite mean energy N — {{ipQ\a^a + 
b^b\'^o))- In order to find physically realizable states, one must impose a con- 
straint on N too, and the Lagrange function becomes 

£ = (7 - fx{{^Po\a^a + b^bl^Po)) - X{{i^o\^o)) , (16) 

II being the Lagrange multiplier for the mean energy. It is useful, in order to 
calculate the solution of this equation, to write the generic state |V'o)) in the 
following way 

oo 

\i^o)) = J2^dJ2^n,d\f^,d)). (17) 
deZ n=0 

The coefficients Cn,d determine the normalized projection of |^)) into the 
eigenspace Hd, whereas the ipdS are the coefficients which combine those pro- 
jections. 

Using Eq. (17), the Lagrange function (16) explicitly shows terms accounting 
for the normahzation of the vectors J2'^=o Cn,d\n, d)) in each Hd, with i/^'') de- 
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noting Lagrange multipliers for the normalization of projections, and rewrites 
as 



d€Z [ n=0 \ n=0 / \ m=0 ) 



oo \ / oo 



\ n=0 / \ m=0 / 

oo oo oo '^ 

X |c„,,p(2n + - J2 Kd\' - ^^'^ E |cn,dP .(18) 

n=0 n=0 n=0 ) 

By taking derivatives of the Lagrange function with respect to c*^ and ip^ 
with the constraints 

oo 

E icn,.r = 1 E iv^'^r = 1 

n=0 deZ 

oo 

YT.\H'\cn4\'i2n+\d\) = N, (19) 

dgZ n=0 

and by rephasing the |(i))'s we arrive at the system 

(2 - X)^Pd - ^d-i - ^d+i - /iE^=o(2n + \d\)\cn4?^d = ^^^^ 
(2 - A)^J - 2(^,V^rf_i + V^d+i^d) - ^^i'l{2n + |rf|) - z/W] c„,rf = 

The second equation in (20) implies that for a fixed d only one coefficient 
can be different from zero, say for the value n, and in this case \cn4\ = 1 Q- 

The first equation of (20) can therefore be rewritten as 

(2 - A)^rf - ^a-i - ^d+i - KMd) + MD^d = 0, (21) 

which allows us to obtain from the second one 

ly^'^ = -Mi^d+i+^d-i)- (22) 

The solutions of Eq. (22) give local minima for the average cost C, and one 
should solve the equation (21) for arbitrary choices of n{d), looking for the 



Indeed, let us suppose Cn^d 7^ foi^ two values ni — tti and 77,2 — p- Then we 
must have (2 — X)tp'^ — 2(7/;rf7/'d-i + V'd+iV'd) — ^V'd(2p + |(i|) — u^'^^ = which imphes 
2fi{m — p)iljd = 0. Since the case = is not interesting and ^pd = would imply 
that the choice of the Cn^d is completely arbitrary and irrelevant, the only possibility 
IS m = p. 
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optimal one. In the case n{d) — we have 



^^^^^V^. = + , (23) 

where a' = —a and A' = Eq. 23 is the recursion equation for Bessel 
functions, with solution given by 

V^, = 7V-V2(^',A')JA'+|d|(^), (24) 



with A/'(//', A') = J2dGZ J\'+\d\ {f/J and with the boundary conditions J^,(2//x') = 
0, i.e. Jy+i = Jy-i- Finally, to obtain the optimal state po one has to rotate 
Ro — |V'o))((V'o| by the unitary transformation (6). 




1 10 100 

N 



Fig. 1. Average cost C as a function of the energy N for the optimal states (dashed 
line). Solid line is the function 1/N'^. The points above the solid line are other 
solutions of equation (21) corresponding to local minima of the cost, distributed 
over lines characterized by a fixed value of the Lagrange parameter (a single 
dotted line is plotted, connecting points for different /x and increasing N). 

In order to obtain the behaviour of the average cost versus the energy we 
numerically solved Eq. (21) with n{d) — 0. This problem can be rewritten as 
the eigenvalue problem Aip — Xip for the matrix A with elements given by 

(A)^,„ = (2 - IJ'\m\)Sm,n - Sm,n+1 " (^m,n-l ■ (25) 



Numerical diagonalization gives the power-law (7 ~ in the range < 
< 1000 with 7 ~ 0.1. This behaviour is plotted in Fig. 1, where also other 
solutions of equation (21) corresponding to local minima of the average cost 
are shown. Since the phase distribution of the optimal state is singly peaked 



7 



we may also write Sip ~ VC ~ j^^^/N, which means that the optimal states 
derived here are at the so-called Heisenberg limit of phase variance. 



3 Conclusions 



In conclusion, wc dealt with the problem of estimating the coupling constant 
of a bilinear interaction. In our approach the coupling constant, which appears 
in the exponent of the time evolution operator, has been treated as a phase 
parameter. The optimal POVM has been derived according to the theory of 
quantum phase estimation [4-6,8]. As noticed in Ref. [2] this resorts to an 
SU(2) estimation problem with the Schwinger two-mode boson realization. 
However, the representation is not irreducible, and a more complicated prob- 
lem is faced. In this sense our results generalize those of Ref. [2], where the 
estimation of SU (2) phase-shifts has been analyzed in irreducible subspaces: 
indeed we found an improved scaling of phase variance versus the total en- 
ergy. The degeneracy of the spectrum of the the Hamiltonian, can be treated 
using the technique of Ref. [6] , and, in this way, the problem is reduced to a 
nondegenerate one with spectrum Z for the phase-shift operator. 



It is important to remark that the optimal POVM depends on the preparation 
state. This is true for every kind of phase estimation problem, but in the 
presence of degeneracy the dependence is crucial. In fact, one has to define the 
optimal POVM as a block diagonal operator, where the invariant subspaces are 
spanned by projections of the input state into the eigenspaces of the generator. 
Prom a practical point of view this means that optimal estimation of the phase- 
shift imposed to a state, needs a measuring device which is adapted to the 
shifted state. We then optimized the input state by minimizing the average 
cost for fixed input energy and found a power law Sip ~ ^/N in a range 



Notice that the law A^~^ is the same as in the optimal phase estimation with 
only one mode [9], however with a much smaller constant 7 (7 ~ 0.1 instead of 
7 ~ 1.36). We think that this phenomenon of improvement of phase sensitivity 
by increasing the number of modes is the same considered in Ref. [10], where 
an exponential improvement versus the number of modes has been estimated 
when increasing the number of modes and the the number of photons per 
mode, jointly in the same proportion. 
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